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Abstract 

This is a continuation of the work on the spectral problem of 
Harper operators. We continue to discuss the local monodromy repre- 
sentation of the Fermi curve using the Picard-Lefschetz formula. The 
density of states over approximating components of the Fermi curve 
satisfies a Picard-Fuchs equation. A mirror map is derived based on 
the Picard-Fuchs equation, which gives a g-expansion of the energy 
level for the electronic model. 



1 Introduction 

In [llj we studied the spectral theory of Harper operators with irrational 
parameters as a discretized model of magnetic Laplacian. On an algebro- 
geometric level, the compactification of its associated Bloch variety, which 
describes the complex energy- crystal momentum dispersion relation, is an 
ind-pro-variety with properties analogous to a totally disconnected space. 
This corresponds to the classical result that the Harper operator's spectrum 
is a Cantor set of zero Lebesgue measure. Further, we expressed the density 
of states of the electronic system as a period integral over Fermi curves, 
which we explicitly computed as elliptic integrals and in terms of the Landen 
transformation. 

In this paper, we continue the investigation of the algebro-geometric prop- 
erties of Harper operators, by focusing on the monodromy of the Fermi 
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curves. The local monodromy can be computed using the Picard-Lefschetz 
formula, while the global monodromy was already discussed in [5] . The den- 
sity of states also has a local monodromy, which leads us to consider the 
ratio of two linearly independent solutions of the Picard-Fuchs equation. We 
prove that the energy level satisfies a Schwarzian equation with respect to the 
complex structure modulus, which is very similar to the form of the mirror 
map for elliptic curves. We also derive a g-expansion for the energy level. 



1.1 Harper operator 

In order to fix the notations, we briefly recall the discretized spectral problem 
of the Harper operator. As a discretized magnetic Laplacian, the Harper 
operator H acting on £^(Z^) is defined as 

H^{m,n)= e-2™"^(m + + e2™"V(^ - 1,^) + 
g-2^i/3m^(^^^ + 1) + e2^*^"'^(m, n - 1) 

where the two unitary operators 

[/V^(m, n) := e-'^^'''''^{m + 1, n) and V^{m, n) := e-2^^'^"7/'(m, n + 1) 

are the so-called magnetic translation operators with phases a and /3 respec- 
tively, so that the Harper operator is of the form H = U + U* + V + V* . 
We are mostly interested in the cases when a and /3 are irrational real num- 
bers. In the theory of noncommutative geometry, the magnetic translation 
operators U, V give rise to a noncommutative torus Ag with 6 = a — [3. 

We consider the spectral problem of the Harper operator in the algebro- 
geometric setting, for distinct primes a and 6, define the Bloch variety asso- 
ciated to the Harper operator ([T]) as 

B:= {(6,6,A)GexC*xC|i/7/; = AV^, 

'?/'(m -|- a, n) = ^i'?/'(m, n), '?/'(m, n + 6) = ^2'^(^, '^)} 

Due to the irrationality of the magnetic fluxes a and /3, the locus determined 
by ([2]) consist of a countable collection of algebraic varieties. More precisely, 
the Bloch variety in this case is an inductive limit of finite dimensional alge- 
braic varieties, we call it the Bloch ind-variety. 
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Lemma 1.1. [77]/ The Block variety defined by the spectral problem of 
the Harper operator ([T]) is an affine ind-variety that can be written as B = 
\Jki&Bk/, where 

Bk,i = {(6, 6, A) eC*xC* X C|P,,,(ei,6,A) = 0}, (3) 

for polynomials 

PM(ei,6,A) = det(M(^'^)-A/), (4) 
with the ab by ab matrix M^'''^^ having entries 
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(5) 



After the change of coordinates by Fourier transform, the Bloch ind- 
variety can be rewritten as 5 = IJ^ ^ ^ ^ Bk/,m,n with nonsingular subvarieties 
as components: 

BuAra,n := {(6,6, A) G C* X C* X C | iV^l(6,6) - A = 0} (6) 
where (fc, £, m, G Z x Z x x Z;, and 

Arfc/ (c (: \ 2iTia{n+ih) c , —2nia(n+£b) f^—l , 2nil3{m+ka) /: , —2TTij3{m+ka) f^—l 

(7) 

Hence it is convenient to work on the approximating components of the 
Bloch ind-variety. 

1.2 Mirror map for elliptic curves 

In this subsection, we review the definition of mirror maps for eUiptic curves. 
In particular, we explain the concepts of Picard-Fuchs equation and the re- 
lated Schwarzian equation by the j-invariant, whose mirror map is the power 
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series expansion of j in terms of q = e^'^*'^, r G HI/S'L(2,Z). More informa- 
tion about mirror symmetry and algebraic geometry, especially its relation 
with Picard-Fuchs equations can be found in |3], [9], [T2] . 

Let X be a Calabi-Yau three-fold and X* its mirror image, the general 
setting is that the B-model correlation function Krrr over X* is depending 
on the complex structure modulus r of X* and such dependence is encoded 
in a Picard-Fuchs equation. In general, a differential equation is called of 
Picard-Fuchs type if all points on are regular or a regular singularity and 
its monodromy group is contained in SL{n,Q), where n is the order of the 
equation and Q is the algebraic closure of Q. 

On the other hand, the Yukawa coupling over X, i.e. the A-model corre- 
lation function, Kjjj has a formal expansion: 



^jjj 



X 



J A J A J + Srf^^dl (8) 
1 — q'^ 



where na is the counting number of rational curves of degree d on X, or the 
instanton numbers for physicists. If t is the complexified Kahler class of X, 
whose imaginary part cubed (Imt)^ is proportional to the volume of X, then 
q = e^'^** as in the above expansion. 

In superstring theory, mirror symmetry simply states that Krrr = Kjjj. 
In order to identify the mirror symmetry, one has to determine the mirror 
map between local coordinates of the complex structure moduli space and 
the Kahler structure moduli space. In particular, the mirror symmetry for 
elliptic curves interchanges the complex structure modulus r of X* with the 
Kahler structure modulus t of X, so that the g-expansion in the B-model (i.e. 
a complex torus) is the same as a formal expansion in terms of q = e^'^*'^. 

Recall that the Schwarzian derivative of a function f{z) is defined by 

Let fi{z) and f2{z) be two linearly independent holomorphic solutions of 



+ Q{z)f{z) = 0. (10) 



Then the ratio g{z) = fi{z)/f2{z) satisfies {Sg){z) = 2Q{z) over the do- 
main on which fi{z) and f2{z) are defined. The Schwarzian derivative has a 
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inversion formula, 



Let us consider a Weierstrass family of elliptic curves = Ax^ — g2X — gs 
embedded in P^, where j = (^1/ ((?| — 27(y'|) is the j-invariant. The correspond- 
ing Picard-Fuchs equation is given by 



(i^TT 1 dn 



+ + 



31j 



-TT = 0. 



df Jdj 144j2(l-j)2 

By a change of variable tt = uv, it is easy to get its Q-form: 



[12) 



d'^u 



+ 



2 



23 



df [16(1 -jT 9f 144j(l-J 



u = 0. 



(13) 



Hence the ratio of two period solutions r = 112/711 = U2/U1 satisfies the 
Schwarzian equation 



4 

+ TT^ + 



23 



9j2 72j(1-j: 



(14) 



Change the variable ^ = 1/j in the Picard-Fuchs equation (IT^ to look 
into its behavior around j = 00, rewrite it as 



A Idy 3ie-4e^ ^ 

Denote the holomorphic functions near ^ = (i.e. j = 00) by 

3ie - 



(15) 



p 



Q 



144(1-0^' 
then the above equation is of the form 

de^idCe^ ' 

meaning the Fuchsian equation (fT2|) has a regular singularity aX j = 00. 
One computes the indicial equation, which is 



r(r- 1) + P(0)r + Q(0) 
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it has a double root at r = 0. Then the two dimensional space of solution is 
spanned by yi{^), which is holomorphic near ^ = 0, and ?/2(0 = Vi^^ ^ + P(0 
for some holomorphic function p{C,)- 

Thus the solution to the Schwarzian equation f|T4|) has the form 

in other words, the corresponding q has the form 

g(0 = e* = e ■ eP'y^ (17) 

This defines an invertible analytic map from a disk |^| < S" to some disk 
|g| < -R sending to 0. 

One can apply the Frobenius method to find out ?/i,y2 and then get the 
g-expansion 

^(q) = J_ = 1 _ 744g2 + 356652g3 - 140361 152g^ + ■ • • (18) 

Then j = j{q) (or its inverse) is the mirror map and its g-expansion is the 
well-known power series by the monstrous moonshine: 

j(q) = 1 + 744 + I96884g + 21493760g^ + 864299970g^ + • • • 
Q 

2 Monodromy of Fermi Curves 

In this section, we recall the definition of the algebraic Fermi curve and dis- 
cuss its monodromy. The local monodromy representation will be derived 
from the Picard-Lefschetz formula and the global monodromy group is al- 
ready observed in [3]. 

We consider the projection vr : — t- C; (^i,^2)A) A, recall that affine 
Fermi curves are defined by Fx{C) := 7r~^(A), then each Fermi curve Fx = 
U/c £ m n F^'^'"^'" is an ind- variety with components 

i^A''''"'^ = {(6,6)Ki(ei,e2) = A} (19) 



has two components 



For each 4-tuple {k,i,m,n), there exist singular fibers at A = 0,4,-4. 

its 

{(ei,e2)|e2-^"(^+^^)6 + e2-^(-+^'^)6 = 0} 
{(6, 6)|e^"*"("+*)6 + e-2-^/3(m+fca)^-^ = 0} ^ ^ 
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pk,e,m,n gjng^lg^r Qj^ly ^g-2™(n+^6)^g-2«/3(m+fca))^ similarly F^;/''"''^ is 

singular only at (_e-27riQ:(n+ffe)^ _g-27rj^(m+A;a)^ g^j^^ they are irreducible curves. 
In addition, these singularities are all ordinary double points. 

As in Lemma 5.1 of |5], the projective closure of each Fermi curve com- 
ponent _p^'^'"*'" in X P^, call it p^'^'"^'"'^ ig elliptic curve for a generic 
A, and the complement p^'^'"^'"^ y pk/,m,n ^ divisor of type (2,2). Then 
we have an elliptic fibration of the projective closure of each Bk/,m,n, call it 

Bk.£,m,ni 

n : Bk,e,m,n ^ F\ 7r-i(A) = (21) 

By the same statement as in [5j, the above fibration gives rise to a stable 
family of elliptic curves with four exceptional fibers, at A = ±4 type /i, at 
A = type I2 and at A = cxd type /g, whose global monodromy group is 
ro(8) n r°(4). The details can be found in ^. 

2.1 Local monodromy and Picard— Lefschetz transfor- 
mations 

The local monodromy of the family of elliptic curves has been discussed in 
Proposition 8.1 of [5]. Here we have a similar result for the family of Fermi 
curve components. 

Let 5 be a vanishing cycle, the Picard-Lefschetz formula is given by 

T(x) = x-{x-6)6 (22) 

for an arbitrary homology cycle x and x-6 is the intersection number between 
cycles with chosen orientation, cf. |13j . 

Lemma 2.1. The local monodromies around A = 4,0,-4 are given by the 
Picard-Lefschetz transformations T4,Tq,T_4, 

T4(7)=7 + 5i, r_4(7)=7-52, To(7) =7-^1 + ^2, (23) 
where 61,62 are vanishing cycles and 7 is an arbitrary homology cycle. 

Proof. First change the variables by ^ = ^'^'^^o'i^+^f') and r] = e27rj/3{m+fca)^^_ 
Then the Fermi curve components become 

(24) 
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Rewrite the above curve as ^rf + (^^ — + l)rj + ^ = 0, whose discriminant 
is A = (^^ — + 1)^ — 4^^. Set A = 0, there exist four branch points 

When A tends to 4, shrink to be an ordinary double point 

at ^ = 1, which corresponds to = (e-2™("+^''), e-2^^^('"+'=")), the cycle 

around jg ^ vanishing cycle, call it 5i. Assume 7 is a cycle such 

that 5i "7 = 1, we rotate 5i once and interchange the points i^-'^)+'^^^ z^ ^ — 
(A-2)-VA^ --4A^ Picard-Lefschetz transformation T4(7) = 7 + 5i. 

When A tends to —4, shrink to be an ordinary double point 

at ^ = -1, which corresponds to (6,6) = (-e-2™("+^^), -e-2^^^('^+'=")), 
the cycle around jg ^ vanishing cycle, call it ^2- 7 is the 

same cycle as above, so ^2 ■ 7 = — 17 we rotate 82 once and interchange 
the points ('^+^)+^'^ +fA ^ — y (a+2)~va +4A^ ^j^^ Picard-Lefschetz transfor- 
mation T_4(7) = 7 — 52- 



When A tends to 0, ('^+^)='=^^^ ±jA ghrink to be one ordinary double point at 
^ = 1 corresponding to (6, 6) = (e~^™'^"'+^^\ _e-27ri/3(m.+A:a)-j g^j^^ ^y^q same 
time (''^~2)='=^^! i:fA ghrink to be the other ordinary double point at ^ = — 1 
corresponding to (6,6) = (-e-2™("+^''), e-2"^(™+'='')). The cycles 5i,52,7 
are the same as above, so by the Picard-Lefschetz transformation Tq{^) = 
7 - 5i + 52- □ 



3 Picard-Fuchs equation 

One of the main results in lHj was that we identified the density of states 
on each component of the Fermi curve with a period integral independent of 
the labeling. In this section, we will give the Picard-Fuchs equation of the 
density of states, and use the Frobenius method to derive a q-expansion of 
the energy level. 

Recall that the density of states (DOS) on each component of the Fermi 
curve p^'^'"^'^ ig of the form 

1 1 9\/h 

where k = |^ is the elliptic modulus. To avoid confusion with the elliptic 
modular lambda function, we write the energy level as = |A| instead of A 
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as in the previous section. The right hand side of the above equation is the 
ascending Landen transformation, which changes the period from r to r/2. 

The complete elhptic integral of the first kind K{k) satisfies a second 
order differential equation 

k{l - k')^ + (1 - 3k')^ -kK = (26) 
ak^ dk 

with regular singularities aXk = —1, 0, 1, oo. In fact, if we change the variable 
k by the elliptic modular lambda function A(r) = k"^ in K, then the above 
differential equation is equivalent to 

which is the Picard-Fuchs equation of the Legendre family = x{x—l){x—X) 
of elliptic curves. 
If we define 

D{k) = 27r^ab DOS = (1 + k)K{k) (28) 
then we have the following result. 

Proposition 3.1. The density of states, or equivalently D{k), satisfies a 
second order differential equation as follows 

d'^D dD 
k{l -k){l + kf-— + (1 - 2A; - + k)— + {k- 1)D = 0. (29) 
dfc^ dk 

We call (12 9 p the Picard-Fuchs equation of the density of states. 

Proof. Recall that the relation between K{k) and the complete elliptic inte- 
gral of the second kind E{k) is the following. 



dK _ E{k) _ K{k) dE _ E{k) - K{k) 
~dk~k{l- F) lik~ k 

We substitute and compute the first derivative of D{k): 



(30) 



In other words, 

k{l-k)'^ = E{k) + {k-l)K{k) (32) 
dfx, 
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Then 

d , ,sdD^ dE .dK , . 



,,d^D , ,,dD E-K E K , E 

- =—+<*- ^Xmi^is) - t) + ^ = m 



Here we use (132!) again to cancel E{k), so 



(34) 



□ 



3.1 Picard— Fuchs equation and Picard— Lefschetz for- 
mula 

Instead of ^, some authors prefer to use Q = k-^. With this notation, we 
can reformulate the Picard-Fuchs equation for the density of states in the 
following way. 

Corollary 3.2. The Picard-Fuchs equation gives rise to a differential oper- 
ator in terms of 9 as 

02 + - — ^ (36) 

Since the monodromy representation for solutions of the Picard-Fuchs 
equation is the same as the geometric monodromy representation, we consider 
the local monodromy of D{k) = K(2\/k/(l + A;)), i.e. the monodromy of 

The corresponding canonical holomorphic form is 

dt dt ds , , 

bj = — = , = , (38) 



where s = — 1/t 



2 _ 4fc 
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Lemma 3.3. The local monodromies around k = 0,1 are given by the Picard- 
Lefschetz transformations Sq, Si, 

^o(7)=7 + 5i, 5i(7)= 7 + 2^2 -2^3, (39) 

Proof. If k approaches 0, then ±2a/A;/(1 + k) shrink to 0, which gives an 
ordinary double point at s = 0. So the cycle around ±2\/k/{l + /c) is a 
vanishing cycle, call it 6i. Let 7 be a cycle intersects with Si so that 5i -7 = 1. 
Hence by the Picard-Lefschetz formula 6*0(7) = 1 + ^i- 

If k approaches 1, then ±2y/k/{l + k) go to ±1, which gives two ordinary 
double points at s = ±1. So the cycles around 2\/k/ (l+/c), 1 and —2\/k/ (1 + 
/c), —1 are vanishing cycles, call them 82, ^3, let 7 be a cycle intersect with 82 
and 8^ such that ^2 ■ 7 = 1 and 53-7 = —1. Similarly, by the Picard-Lefschetz 
formula 5*1(7) = 7 + 2(^2 - 253. □ 

By the relation k = if — )■ 4, or equivalently A — )■ ±4, then 

A; — )■ 0, and this gives an ordinary double point as in section [2]; if -E — ?■ 0, 
or equivalently A — )■ 0, then A; — )■ 1, and this also gives rise to two ordinary 
double points as before. 

Notice that, in the elliptic curve model of two dimensional Ising model, 
the point k = 1 corresponds to the critical temperature of the Ising model. 



3.2 The Q-form and the Schwarzian equation 

Monodromy of the Picard-Fuchs equation of density of states means that if 
we go around some singularity then the period solutions (density of states) 
will change according to the monodromy representation, which forces us to 
consider the ratio of the period solutions modulo the monodromy group as 
in the theory of elliptic curves. In the following we change the Picard-Fuchs 
equation f l29p into its Q-form and consider the related Schwarzian equation. 

Lemma 3.4. By change of variable D = UV , the Picard-Fuchs equation is 
equivalent to 

^ + Ji±^^^0 (40) 
the corresponding differential operator is 



2 



+ (41) 
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Proof. Based on (12^ . first define V such that 
and the second derivative is 

7P + 2^ _^ (43) 

Then by substituting into the Picard-Fuchs equation and canceling out V 
and the middle term, it is easy to get the Q-from, where 

□ 

Proposition 3.5. The energy level E satisfies a Schwarzian equation: 

Proof. Let Di{k), D2{k) be two linearly independent period solutions of the 
Picard-Fuchs equation fl29l) . and let Ui{k),U2{k) be the corresponding solu- 
tions of the Q-from fHOj) . define r = D2/D1 = U2/U1, then r satisfies the 
Schwarzian equation 

With D{k) = K{2y/k/l + k), the complementary modulus is given by the 
descending Landen transformation 



k' 



Ak 1-k 



(l + fc)2 1 + k 
then we can choose r (i = is added for convenience) as 



r = ^Ki^)/Ki^) = - (47) 
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where r G H/F is the complex structure modulus of elliptic curves. Then 
the above Schwarzian equation is expressed as 



So by the inverse formula, 



,2r - A;'2(l + P)2 



{k, r} = -k'irYir, k} = ' (49) 

On the left hand side, we replace k = {A — E) / {A + E) hy E since the 
Schwarzian derivative is GL{2, M) invariant, on the right hand side, we change 
the variable by A(r) = k"^. 

where q = e™"^ is the nome. □ 
3.3 Energy level and expansion 

In this subsection, we derive a g-expansion for the energy level based on the 
Frobenius method on Picard-Fuchs equations. 

The Picard-Fuchs equation fl29]) is equivalent to 

d^D l l-2k-e dD 1 -fe 

l^^k l-P 'dk^¥{l + ky ~ ' *~ ^ 

which has a regular singularity at A; = 0. We want to find out two linearly 
independent solutions around 0. If we set 

then the indicial equation is given by 

r(r - 1) + rp(0) + g(0) = = 0. (53) 
It has repeated roots r = 0, which means the local solutions are given by 

oo oo 

D-^ = J2dnk'', D2 = Dilnk + ^Cnk'' (54) 

n=0 n=l 

by the Fuchs's theorem. 
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Lemma 3.6. Under the normalization -Di(O) = 1, the local solutions are 

(55) 



Di — 1 + k + \k'^ + \k'^ + ^k^ + ^k^ + ^k^ + ^k'^ + • • • 

D2 — Diln k -\- -jk"^ -\- -jk"^ + Tqs^^ Tor^^ ~^ T^?m^^ ~^ Tvm^'^ ~^ 



4'" ' 4'" ' 128 ' 128 ' 1536 ' 1536 

Proof. 

oo oo 
n=l n=2 

Plug into the Picard-Fuchs equation, we get 

En=0 - + En=l ndn{-k^+^ - 3k^+' - k^ + k^-') + 

En=2 ^(^ - l)(^n(-A;'*+^ - + A;" + A;"-^) 
It is equivalent to 

c?i — (io + (c^o ~ di)k — Sdik"^ — dik^+ 

En=i -(^ + 3)(n + l)dn+ik^+^ - (n + l)2rf„+iA;-+=^ 
From this, we have 

di — do, d^ — \do, dz — |rfo, 
{n + l)^d„+i = [1 + 2n - n'^]dn + (n^ - 2)dn_i + (n - 2)2d„_2 



• rfo (56) 



It turns out that 

d2i - d2i+i - [ (20!! 

If we apply the normalization so that -Di(O) = 1, then do — 1. 
Then we look at D2 — Diln k + En=i ^n^", 

D'^ = D[lnk + iDi + En=i ncn k^-\ 

D'i = D'iln k + \d[ - ^Di + E„=2 n{n - l)c„ k^'^ 

Plug into the Picard-Fuchs equation, 

k{l-k){l + kf [D'iln k + lD[ - ^D, + Y.n=2 <^ - l)cn 
(1 - 2A; - A;2)(l + A;)[D;/nA; + + E„=i + 
(A;-l)[Di/nA; + En=iCnA;"] 
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Cancel the terms with Ink since Di satisfies the Picard-Fuchs equation, the 
other terms can be used to determine €„. 



2(1 -k){i + kfD[ + ik-m+k? j)^ ^ (i-2fc-fc^)(i+fc) ^^^ 

{k - 1) En=i Cn A;" + (1 - 2A; - k''){l + k) nc„ A;"-i+ 
A;(l - A;) (1 + A;)^ En=2 ^(^ - l)cn A;"-^ 

Use the power series of Di , this can be reduced to 

En=o(2(^ + ^)dn+i - 2dn + {u + l)c„+i)A;"+ 
(2(n + l)c/„+i - 2dn + {n + l){n + 2)c„+2 - (n + 2)c„+i)A;"+i+ 
(-2(n + l)rf„+i + (n + l)(n + 2)c„+2 - (3n + 2)c„+i)A;"+2+ 
(-2(n + l)dn+i - {n + l){n + 2)c„+2 - (n + l)c„+i)A;"+3^ 
-(n + l)(n + 2)c„+2A;"+4 



Then we have 



Ci = 0, C2 = 2, C3 - ^ 



4 ) "-'o 4 ) 

(n + l)2cn+i = (1 + 2n - n2)c„ + (n^ - 2)c„_i + (n - 2)2cn-2 

-2(n + l)dn+i - 2(n - l)rf„ + 2ndn-i + 2{n - 2)dn-2 



□ 



If we do the long division, then 

7^' + ^^' + ^^' + --- (57) 



12n=2 ^nk"" _ 1 , 2 , 13 4 23 g 



En=odnk- 4 128 384 
Exponentiate the above quotient, we have 

exp{J2 CnkVYl dnkn = 1 + + + + • • • (58) 

n=2 n=0 

Theorem 3.7. T/ie energy level has a q- expansion 

39 11 

£; = 4 - 8? + 8?^ - + 4g^ + —q^ (59) 

16 8 

where q — e^^^^^. 
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Proof. Set t — D2/D1, we have the local coordinate q — as a, power series 
q^k- exp{J2 Cnk^/Y. ^k + ^e + ^k' + ^k' + ■■■ (60) 

n=2 n=Q 

This is the mirror map q — q{k), its inverse is given by 

k = q-\q' + Y^g' + • • • (61) 
Since the elhptic modulus k is defined by A; = it is easy to see that 

1+^= i+q-l'l^+Tkl'+- (62) 

= 4(1 -2? + 2g^- 1^3 + ^4 _ 1^5+ 11^6 + ...) 

□ 

Thus, we see that the behavior of the map E{r) and its g-expansion E{q) 
discussed above have a close similarity to the properties of the mirror map. 
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